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ABSTRACT: CRBP1 (cellular retinol binding protein 1) and CRBP3 (cellular retinol binding protein 3), are important components of 
the retinoid signaling pathway and take part in vitamin A absorption, transport and metabolism. Based on the role of vitamin A in 
chicken laying performance, we investigated the polymorphism of CRBP1 and CRBP3 genes in 349 chickens using single strand 
conformation polymorphism and DNA sequencing methods. Only one polymorphism was identified in the third intron of CRBP1, two 
polymorphisms were detected in CRBP3; they were located in the second intron and the third intron respectively. The association studies 
between these three SNPs and laying performance traits were performed in Erlang mountainous chicken. Notably, the SNP g. 14604G>T 
of CRBP1 was shown to be significantly associated with body weight at first egg (BWFE), age at first egg (AFE), weight at first egg 
(WFE) and total number of eggs with 300 age (EN). The CRBP3 polymorphism g.934C>G was associated with AFE, and the 
g.l324A>G was associated with AFE and BWFE, but none of these polymorphisms were associated with egg quality traits. Haplotype 
combinations constructed on these two SNPs of CRBP3 gene were associated with BWFE and AFE. In particular, diplotype H2H2 had 
positive effect on AFE, BWFE, EN, and average egg-laying interval. We herein describe for the first time basic research on the 
polymorphism of chicken CRBP1 and CRBP3 genes that is predictive of genetic potential for laying performance in chicken. (Key 
Words: CRBP1, CRBP3, Laying Performance, Egg Quality, Gene Polymorphism, Erlang Mountainous Chicken) 



INTRODUCTION 

Retinoids, metabolites of vitamin A (retinol), play some 
essential physiological roles in both vertebrates and 
invertebrates through control of cell differentiation, 
proliferation and apoptosis (Linan-Cabello et al., 2002; 
Reichrath et al., 2007). Besides these functions, retinoids 
also affect many other vital processes including vision 
(Rando et al., 1994), reproduction (Eskild et al, 1994), 
immune function (Ross et al., 1994), growth and 
development (Hofmann et al., 1994) and hematopoiesis. 
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Since animals cannot synthesize retinoids, vitamin A is 
usually obtained from the diet (Novak et al., 2008; Harrison 
et al., 2012). However, once the vitamin A deficiency 
occurs, many biological functions within the body will be 
affected or destroyed. For example, in chickens, lack of 
vitamin A will affect the egg production, egg quality, 
production and immune response (Lin et al., 2002; Abdalla 
et al., 2009; Xiao et al., 2011). In rats, the case of vitamin A 
deficiency has been reported to lead to congenital 
malformation in embryos (Warkany et al., 1948; Li et al., 
2010). 

The hydrophobic molecular structure of vitamin A and 
its biological derivatives (the retinoids) means that they 
must be chaperoned both intracellularly and extracellularly 
by retinoid binding proteins (Newcomer, 1995). Serum 
retinol-binding protein and cellular retinol-binding protein 
(CRBPs) respectively is the extracellular and intracellular 
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retinol binding protein (Guo et al., 2004). At present, four 
types of CRBPs (CRBP1, CRBP2, CRBP3, and CRBP4) 
have been reported. They are highly conserved among 
mammals (Ong, 1994; Zhang et al., 2012), but have distinct 
tissue distributions. CRBP1 and CRBP4 are widely 
distributed throughout the body; meanwhile, CRBP1 is 
expressed at high levels in the liver, testes, eye and kidney 
(Piantedosi et al., 2005). CRBP2 is expressed only in the 
small intestine (Yamaguchi et al., 2009), whereas CRBP3 is 
restricted to the heart, adipose tissue and skeletal muscle 
(Caprioli et al., 2004). Of these four, CRBP1 is thought to 
be essential for vitamin A metabolism and synthesis of 
retinoic acid (Kuppumbatti et al., 2000; Noy, 2000). 
Although these members of CRBPs have a different 
expression features, they all affect the activity of enzymes 
involved in retinoid metabolism (Napoli, 1993). In past 
studies, some researchers found that excess or deficiency in 
retinol were associated with egg laying, egg size, embryo 
survival and hatchability (Thompson, 1969; Boily et al., 
1994). In addition, Xiao et al. (2011) and Yin et al. (2013) 
reported that the variation of CRBP2 and CRBP4 had 
significant effects on egg production traits. Therefore, the 
intracellular retinol binding protein, CRBP1 and CRBP3 
may play important role in transport, metabolism and 
homeostasis of retinol. This makes CRBP1 and CRBP3 as 
good candidate genes for enhancing reproductive traits in 
chickens. Therefore, to prove this hypothesis, the objective 
of the current study was to characterize polymorphisms in 
the partial exons and introns of the CRBP1 and CRBP3 
gene, and determine their association with laying 
performance traits in Erlang mountainous chickens. 

MATERIAL AND METHODS 

Animals and DNA preparation 

Two lines of SD02 (Shandi chicken) (n = 188) and 
SD03 (n = 161) were enrolled in this study. These chickens 
are same as used in Xiao's research (Xiao et al., 2011). 
SD02 and SD03 are new strains that all have the blood of 
Erlang mountainous chickens (Erlang mountainous chicken, 
an indigenous chicken breed distributed in Tianquan of 
Sichuan province, widely adaptable, and have delicious and 
nutritious meat) and all were developed on an experimental 
farm for poultry breeding at the Sichuan Agricultural 
University (Ya'an, China). Meanwhile, there are differences 
in the specifics of these two strains. SD02 has yellow 
partridge plumage, blue shanks, white skin and heavy body 
weight, while SD03 has normal body weight. During the 
growth period, all birds had free access to food and water 
ad libitum under the same temperature and lighting 
conditions. 

Six traits were measured and recorded: the age at first 
egg (AFE), highest continuous egg days (HCD), body 



weight at first egg (BWFE), weight at first egg (WFE), total 
number of eggs with 300 ages (EN), and average egg-laying 
interval (AEI). Approximate 3 mL blood per chicken was 
obtained from the wing vein and kept in a vacuum tube 
containing anticoagulant acid citrate dextrose (ACD). All 
samples were delivered to the laboratory in an ice box. The 
genomic DNA was extracted using a standard phenol- 
chloroform method. The DNA samples were dissolved in 
TE buffer and were stored at -20 until used. 

PCR amplifications 

The primers were designed on the basis of DNA 
sequence of the CRBP1 (Accession no. NC_006096.2) and 
the CRBP3 (Accession no. NC_006088.2) using the oligo 
nucleotide design tool Primer 5.0 software (Table 1). Primer 
synthesis was completed by Shanghai Yingjun 
Biotechnology Co. Ltd. (Shanghai, China). PCR reactions 
were performed using the Gene Amp PCR System 9700 
(Bio-Rad, Hercules, CA, USA) thermal cycler in a final 
volume of 10 uL reaction containing 10 to 100 ng genomic 
DNA, 5 (iL 2xTaq PCR MasterMix (including Mg 2+ , dNTPs, 
Taq DNA polymerase; Beijing TIAN WEI Biology 
Technique Corporation, Beijing, China), 3.4 |iL ddH 2 0 and 
0.4 (iL of each primer (10 pmol/uL). The cycling protocol 
was 4 min at 94°C, 35 cycles of denaturing at 94°C for 30 s, 
annealing at X°C (Table 1) for 30 s, extending at 72°C for 1 
min, with a final extension at 72°C for 8 min. The PCR 
products were characterized in 1% agarose gel 
electrophoresis. 

Single strand conformation polymorphism and DNA 
sequencing 

Genetic variants in the CRBP1 and CRBP3 genomic 
sequence were analyzed by using the PCR-SSCP (single 
strand conformation polymorphism) method. Each PCR 
product was diluted in denaturing solution (95% formamide 
deionized, 0.05% of bromophinol blue, 0.25% xylene 
cyanole and 10% glycerol) and was denatured at 99°C for 
10 min, then quickly chilled on ice for 5 min. Three 
microliters of mixture was resolved on 12% polyacrylamide 
gel (polyacrylamide:bisacrylamide = 39:1) using lxTBE 
buffer and run for 9 to 14 h at 130 to 150 voltage. The gel 
was washed in 10% ethanol for 10 min, then in deionized 
water for 1 min, followed by an incubation for 10 min in 
0.1% silver nitrate, and ended with an incubation in 1.5% 
sodium hydroxide, 0.01% sodium borohydride, and 0.4% 
(v/v) formaldehyde for 10 min. The bands on the gel were 
visualized by gel imaging system. After the polymorphisms 
were detected, each of the DNA bands on the Single strand 
conformation polymorphism (SSCP) gel was extracted and 
the PCR products of the different electrophoresis patterns 
were sent for sequencing in both directions (repeated three 
times) in an ABI 377 DNA analyzer Applied Biosystems 
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Table 1. Detailed information of the primers used for the single nucleotide polymorphism scanning of CRBP1 gene and CRBP3 gene 



Gene name 



Primer 



Sequences (5 '-3') 



Annealing 
temperature (°C) 



Fragment 
length (bp) 



CRBP1 


Exonl 


F 
R 


ACAATTACCTGAAGTAGGGA 
ACAAAAATGCTTTATGACCAT 


47.6 


189 




Exon2 


F 
R 


AAGCTCTTTTGGTGGTAGATG 
CATGCATTTGCGATCA 


45.2 


234 




Exon3 


F 
R 


ACATGCCCGTTTTCTTAGAC 
AGCACTGCATCTAGCCGTTTC 


58 


185 




Exon4 


F 
R 


ATTATTTTCTCTCCCTAGGAA 
GGGGTTTTATTTGCAATTAT 


52.5 


255 




Intron 1 


F 
R 


ATGGTGATGTACTTGATGCAG 
TGGGAGAAATGAGGTGAAC 


51 


220 




Intron 2 


F 
R 


TGTATTTGTGGCTATTCGTAA 
AGCTTTTGAGCACCTGATTGT 


50 


307 




Intron 3 


F 
R 


TAAATTTATAGTTGCGGTCAT 
AGTTCCAGGTGACGGGTGAG 


50.7 


320 


CRBP3 


Exonl 


F 
R 


GATGGCGTAATGGCTCTAAAG 
CATCCCCGTTGAGACTATCA 


48.5 


104 




Exon2 


F 
R 


AGGTGCTCTGTCTCGGCTTTA 
AATGACCGAAGGGTTGGTTGC 


48 


199 




Exon3 


F 
R 


GAGCCCCCTGTCCAAAGAGTG 
GGCACTTGCGCCCATC 


54 


245 




Exon4 


F 
R 


ACTGCTTTTGTTGGGTCCACG 
ACCAGGCAAGCTGCATGGG 


56.9 


176 




Exon5 


F 
R 


GCCCTTCCTACCCTCGG 
CGCTGCCATCGAGTTTGCTTC 


60.8 


229 




Intron 2 


F 
R 


ACTCCTCCCTCGTTCTT 
CAGGCTCCTCCACTCTT 


52.5 


193 




Intron 3 


F 
R 


GGGATGAACCAGAATG 
CCTTGGGACTGATGCT 


51.1 


213 




Intron 4 


F 
R 


TTGCCTGGTGGTTGG 
GCGGCGGATGTTTAT 


53.2 


137 



and the sequences were determined with DNASTAR 
software (version 7.1, DNASTAR, Inc., Madison, WI, 
USA) and blast in NCBI (National Center for 
Biotechnology Information). 

Egg quality 

Egg quality parameters of egg weight, eggshell 
thickness, eggshell strength and yolk color were measured 
using 150 eggs at 300 day. Eggshell thickness was 
measured with a vernier caliper to the nearest 0.01 mm. Egg 
weight (Wg) was determined with an electronic balance to 
the nearest 0.001 g. The eggshell strength, yolk color and 
albumin height (H) was measured using an egg quality 
measuring tester (EMT 5200, Robotmation Co. Ltd., Tokyo, 
Japan). Then, the egg shape index was calculated as egg 
diameter divided by an egg length. Haugh unit was 
determined using the following equation: 



Statistical analyses 

Genotypic and allelic frequencies and Hardy-Weinberg 
equilibrium were directly calculated. Population genetic 
index such as polymorphism information content (PIC) was 
performed by PopGene software (Version 1.3.1, University 
of Alberta, Edmonton, Canada) (Yeh et al., 1999). Sequence 
analysis was performed with DNASTAR software (version 
7.1). 

Haplotypes of the CRBP3 gene were obtained based on 
the identified SNPs in 349 chicken samples using the 
PHASE 2.0 software (University of Chicago, Chicago, IL, 
USA) (Stephens et al., 2001). 

The association analyses between the single marker 
genotypes and laying performance traits were performed by 
the PROC general llinear model (GLM) procedure of SAS 
6.12 (Statistical Analysis Systems Institute Inc., Cary, NC, 
USA). The model is as follows: 



Haugh unit = 100xlog(H-1.7xWg°' 37 +7.57) 



Y, = n+G 1 +S J +(GS) 1J +e 1J 
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where Yy was the trait measured in animal, |i was the 
population mean, was the effect of the genotype for SNP 
in the CRBP1 and CRBP3 gene, Sj was the effect of sex, 
(GS)ij was the interaction effect of genotype i by sex j, ey 
was the residual random effect associated with animal. Data 
were presented as least square means with associated 
standard error. 

Analyses of the associations between the combined 
haplotype and laying performance traits were carried out. 
The model was similar to that of single marker association 
analysis. 

The additive genetic value at each locus was determined 
using the equation of (BB-AA)/2, where A was the first and 
B was the second allele of the analyzed markers. Then, the 
dominance genotypic value for each locus was equal to the 
difference between the mean of the heterozygotes and the 
midpoint of the two homozygotes (AB-[AA+BB]/2). 

RESULTS 

Genotype and allele frequencies in chicken populations 

In the present study, we first used the PCR-SSCP 
method to check if there were different bands of PCR 
product for each primer, and then applied the direct 
sequencing method to find the possible variations in the 
different bands of the product. In brief, combining different 
bands on the SSCP gels (see Materials and Methods 
section) and sequencing, we detected one SNP in the 
CRBP1 gene and two SNPs in the CRBP3 gene. We named 
these polymorphisms relative to their respective positions 
on sequences of chicken CRBP1 and CRBP3 genes 
(GenBank accession number NC_006096.2 and 
NC_006088.2, respectively). The locations of these SNPs 
are as follows: for CRBP1 gene, one in intron 3 
(g.l4604G>T); for CRBP3 gene, one in intron 2 
(g.934C>G) and one in intron 3 (g.l324A>G). 

The strains effect was not significant; therefore, the data 
from the two varieties of Erlang mountainous chicken were 
pooled and analyzed together. The results of the genotypes 
and alleles of the identified SNPs in the CRBP1 and CRBP3 
genes are shown in Table 2. For CRBP1 gene, the allele G 
was the dominant allele and presented the highest allele 
frequencies in Erlang mountainous chicken populations 
(average = 54.17%). The PIC was 0.3730. The overall 



observed genotypic frequencies deviated from the Hardy- 
Weinberg equilibrium (X 2 = 16.84, p<0.05). For CRBP3 
gene, considering the first polymorphic site, allele g.934G 
was the most frequent in the chickens we analyzed 
(64.72%). The most frequent genotype was g.934GG 
(48.78%), and the least frequent was g.934CC (19.34%). 
The PIC was 0.3521. The overall observed genotypic 
frequencies deviated from Hardy-Weinberg equilibrium (X 2 
= 20.96, p<0.05). Similar with g.934C>T, the second SNP 
showed a higher frequency of the g.l324G allele (62.35%). 
The g.l324AG genotype was the most frequent (48.25%), 
and the g.l324AA genotype was the least (13.53%; PIC = 
0.3587). The overall observed genotypic frequencies were 
in Hardy-Weinberg equilibrium (X 2 = 1.20, p>0.05). 

Association with laying performance traits and egg 
quality traits 

The results of the GLM analysis of association between 
the CRBP1 and CRBP3 gene polymorphisms and laying 
performance traits are summarized in Table 3. As shown in 
Table 3, for CRBP1 gene, the DNA variation at locus g. 
14604G>T SNP revealed a highly significant association 
with AFE, BWFE, WFE and EN (p<0.05). The chicken of 
genotype "TT" had lower AFE, BWFE, WFE, and EN than 
those of genotype "GG" and "GT". 

For CRBP3 gene, the SNP at g.934C>G was 
significantly associated with AFE (p<0.05) (Table 3). The 
chicken of genotype "CC" had the lowest AFE. Meanwhile, 
the polymorphism g.l324A>G showed that the AFE of 
chickens with the AA genotype was significantly lower than 
those with the AG and GG genotypes, respectively (p<0.05). 
The BWFE values of chicken with the AG genotype were 
significantly lower than those with the AA and GG 
genotypes, respectively (p<0.05). No significant difference 
was detected for other laying performance traits. 

However, no significant difference was detected for egg 
quality traits in either SNP of CRBP3 gene; meanwhile, for 
CRBP1 gene, the genotypes were only significantly 
associated (p<0.05) with egg-shape index (data not shown). 

Association of diplotypes with chicken laying 
performance and egg quality traits 

Since there were two SNPs in the CRBP3 gene, the 

potential effect of the haplotypes g.934C>G and 



Table 2. Genotypic and allelic frequencies 


of single nucleotide polymorphisms of CRBP1 and CRBP3 gene 




Gene 


SNP 


Gene frequency 




Genotype frequency 




X 2 


PIC 


G 


T 


GG 


GT 


TT 


CRBP1 


g.l4604G>T 


0.5417 


0.4583 


0.3313 


0.4209 


0.2478 


16.84 (p<0.05) 


0.3730 


CRBP3 


g.934C>G 


C 


G 


CC 


CG 


GG 


20.96 (p<0.05) 


0.3521 






0.3528 


0.6472 


0.1934 


0.3188 


0.4878 








g.l324A>G 


A 


G 


AA 


AG 


GG 


1.20(p>0.05) 


0.3087 






0.3765 


0.6235 


0.1353 


0.4825 


0.3822 







SNP, single nucleotide polymorphism; PIC, polymorphism information content. 
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Table 3. Association of the SNP of CRBP1 and CRBP3 gene and laying performance traits 

Genotype (N) Effect 



Trait 


GG (134) 


GT(164) 


TT(51) 


Additive 


Dominance 


AFE (d) 


176.15 A +10.07 


177.91 A ±10.68 


162.33 B ±8.19 


-6.91 


8.67 


BWFE (g) 


1,897.54 A ±246.24 


1,868.41 A +251.01 


2,039.71 B +236.40 


71.09 


-100.22 


WFE(g) 


40.29 A ±5.93 


41.18 A +6.23 


36.74 B +5.00 


-1.77 


2.67 


EN (egg) 


78.16 A ±17.30 


78.70 A +17.80 


91.67 B +15.84 


6.76 


-6.21 


HCD (d) 


7.35+3.51 


7.89+4.56 


8.43+4.91 


0.54 


0.00 


AEI (d) 


1.70+0.72 


1.69+0.96 


1.46+0.54 


-0.12 


0.11 






CRBP3 




















Trait 


CC (67) 


CG(113) 


GG (169) 






AFE (d) 


169.25 A +11.79 


174.50 Ba +10.80 


177.52 Bb ±10.81 


4.13 


1.12 


BWFE (g) 


1,954.18+262.39 


1,914.96+258.90 


1,878.08+243.17 


-38.05 


-1.17 


WFE(g) 


39.78+6.00 


39.79+5.92 


40.61+6.29 


0.42 


-0.41 


EN (egg) 


84.07+20.13 


81.02+18.01 


78.50+16.73 


-2.78 


-0.26 


HCD (d) 


8.43+5.18 


7.88+4.64 


7.41+3.50 


-0.51 


-0.04 


AEI (d) 


1.64+0.71 


1.69+1.04 


1.66+0.70 


0.01 


0.04 




AA(48) 


AG (168) 


GG(133) 






AFE (d) 


169.15 A +11.88 


175.35 B +12.98 


176.56 B +8.00 


3.71 


2.49 


BWFE (g) 


1,997.29 A +259.87 


l,862.68 Ba +254.99 


l,924.17 b +238.20 


-36.56 


-98.05 


WFE(g) 


39.54+4.72 


39.65 ±6. 11 


41.09+6.49 


0.78 


-0.67 


EN (egg) 


83.35±17.27 


80.28+19.79 


79.44+15.49 


-1.96 


-1.12 


HCD (d) 


7.19+3.62 


7.95+4.58 


7.73+4.03 


0.27 


0.49 


AEI (d) 


1.61+0.72 


1.68+0.86 


1.65+0.81 


0.02 


0.05 



N, number of individuals; AFE, age at first egg; BWFE, the body weight at first egg; WFE, the weight at first egg; EN, the total number of eggs with 300 
age; HCD, the highest continuous egg days; AEI, average egg-laying interval. 

Values are presented by the least squares means+standard error. Values marked with different superscript differ significantly (p<0.05). 



g.l324A>G on laying performance and egg quality traits 
was analyzed. Four possible haplotypes were discerned: 
haplotype 1 = C-G with a frequency of 20.31%, haplotype 
2 = C-A with a frequency of 15.07%, haplotype 3 = G-G 
with a frequency of 41.87%, and haplotype 4 = G-A with a 
frequency of 22.75%. The mixed model analysis indicated 
that there were significant associations between the 
haplotypes and partial laying performance traits (Table 4). 



Haplotypes were associated with BWFE (p<0.05) and 
highly associated with AFE (p<0. 01). There was no 
significant association between the haplotypes and other 
traits. Significant and suggestively dominant effect of the 
diplotype H2H2 was observed for AFE, BWFE, EN, and 
AEI, the diplotype H3H3 had higher WFE, and the 
diplotype H1H2 had higher HCD. Notablely, diplotype 
H3H4 had a negative effect on AFE and BWFE, H1H1 had 



Table 4. Associations of diplotypes of CRBP3 gene with laying performance traits 



Traits 



Diplotypes 


N 


AFE** 


BWFE * 


WFE 


EN 


HCE 


AEI 


H1H1 


23 


174.78+7.85 


1,929.78+258.00 


39.24+4.81 


76.30±15.05 


7.17+3.86 


1.77+0.67 


H1H2 


32 


167.59+13.49 


1,916.56+257.87 


39.45+6.92 


87.44+23.59 


10.03+6.27 


1.64+0.81 


H1H3 


35 


175.20±8.06 


1,960.29+227.90 


40.54+6.85 


81.54+18.41 


7.71 ±3. 54 


1.68+1.20 


H2H2 


12 


163.08+9.14 


2,101.25+252.21 


41.71+5.43 


90.00±14.56 


6.58+2.50 


1.36+0.39 


H2H3 


64 


176.09+11.90 


1,867.03+267.06 


39.50+5.70 


80.14±19.07 


7.97+5.05 


1.72+1.04 


H2H4 


14 


165.50±6.94 


2,020.71+259.12 


39.26+4.40 


83.71+11.63 


7.93±5.34 


1.62+0.48 


H3H3 


75 


177.73+7.93 


1,905.60+237.87 


41.92+6.69 


79.43+14.14 


7.91+4.32 


1.60+0.60 


H3H4 


72 


178.14+12.47 


1,834.86+241.89 


39.88+6.16 


77.22±17.97 


7.01+2.59 


1.66±0.71 


H4H4 


22 


174.77+13.39 


1,925.68±253.67 


38.53+4.32 


79.50+20.78 


7.05+2.80 


1.75+0.94 



Single "*" means that there is significant different between least mean squares for a certain trait (p<0.050); Double "**" mean there is great significant 
different between least mean squares for a certain trait (p<0.01). 

N, number of individuals; AFE, age at first egg; BWFE, the body weight at first egg; WFE, the weight at first egg; EN, the total number of eggs with 300 
age; HCE, the highest continuous egg days; AEI, Average egg-laying interval. 
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the highest reducing effect on EN, and H4H4 had the 
highest reducing effect on WFE. However, there were no 
significant association between the haplotypes and all egg 
quality traits (data not shown). 

DISCUSSION 

Vitamin A (retinol) is a fat-soluble vitamin that has 
received attention with respect to its important role in many 
physiological functions in poultry. Furthermore, the level of 
vitamin A can affect egg production, feed intake, egg 
production rate (Kucuk et al., 2003; Shao et al., 2010) and 
egg quality. However, due to vitamin A and its metabolites 
being hydrophobic, they are bound to CRBPs (type 1, 2, 3 
and 4) in the cytoplasm (Mezaki et al., 2012). Therefore, in 
this study, SNP analysis of two different types of cellular 
retinol-binding protein, CRBP1 and CRBP3, was performed 
to assess whether variation in their sequence could affect 
the egg production and egg quality in Erlang mountainous 
chicken. 

This study, which was designed as the first step in 
detecting genetic markers for egg production and egg 
quality of chicken, uncovered variation in the chicken 
CRBP1 and CRBP3 gene sequence that has not been 
reported in previous research and may provide information 
for a future marker assisted selection (MAS) program. In 
this study, we detected one SNP (g. 14604G>T) in CRBP1 
gene and two SNPs (g.934C>G and g.l324A>G) in CRBP3 
gene. Association analysis revealed that the g. 14604G>T of 
CRBP1 gene was significantly associated with AFE, BWFE, 
WFE, and EN, but not with other traits. Moreover, the 
g.934C>G and g.l324A>G of CRBP3 gene were all 
significantly associated with AFE, meanwhile, the BWFE 
values of g.l324A>G in CRBP3 gene with the AG genotype 
were significantly lower than those with the AA and GG 
genotypes. Several studies addressed the effect of the 
vitamin A on laying rate and egg production. Since CRBP1 
and CRBP3 are members of the CRBPs protein family, the 
mutations might affect the vitamin A metabolism and 
transportation, which may influence cell proliferation, 
differentiation and development. We believe that this is 
reason for the CRBP1 and CRBP3 gene mutations being 
associated with partial egg production traits. 

On the other hand, for egg quality, no significant 
difference was detected for egg qualities (egg-shape index, 
average egg weight, eggshell strength, eggshell thickness, 
haugh unit, yolk color and the percentage of yolk) in 
g.934C>G and g. 1324A>G of the CRBP3 gene in our study. 
This result is in agreement with those reported by Hao et al. 
(2000) and Wang (2008). While, Shao et al. (2010) and Ren 
et al. (2007) reported that vitamin A has a significant 
influence on average egg weight. Moreover, the genotype of 
g. 14604G>T in the CRBP1 gene was associated with egg- 



shape index, the opposite result was found in Wang's report 
(Wang, 2008). In addition, these three SNPs of CRBP1 and 
CRBP3 which were obtained from intron 2 and intron 3 
were synonymous and caused no change in amino acids; the 
mutation still had an effect on partial egg production traits. 
So, we think that there are two possible explanations could 
be offered. One is the mutations could affect CRBP1 and 
CRBP3 function by altering the mRNA stability (Capon et 
al., 2004). The other is that the nucleotide substitution or 
the frame-shift mutation of the genetic mutation may be 
able to change the amino sequence of the target gene, or be 
terminated without producing peptide synthesis of complete 
peptide chains. 

However, some research reported that associations of 
haplotypes with important economic traits were more 
accurate than those of single SNPs (Daly et al., 2001; 
Zhang et al., 2008). This result implied that there was an 
interaction between different SNPs, and that the haplotypes 
generally provided more information than SNPs (Stephens 
et al., 2001). Therefore, haplotype play an important role in 
association studies. We reconstructed haplotypes based on 
two common CRBP3 gene SNPs and discerned an 
association between the diplotypes and egg production traits. 
Particularly, H2H2 diplotype was found to contribute more 
to AFE, BWFE, EN, and AEI than the other diplotypes. 

In general, this study is the first to show an association 
between a DNA variation in the CRBP1 and CRBP3 genes 
and egg production and egg quality in chicken. Associations 
were found for one CRBP1 SNP and two CRBP3 SNPs with 
several egg production traits but they had no effect on egg 
quality. To contribute to genetic improvement of egg 
production and egg quality of chicken, future studies with 
larger sample size and with chicken with phenotypic 
information are necessary to better understand the 
regulatory mechanisms of polymorphism and parent-of- 
origin effect on gene expression of CRBP1 and CRBP3. 
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